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■  A  numerical  method  for  analyzing  large  deflections  of  shells  of  revolution  under  axlsymmetric 
loading  with  arbitrary  creeo  characteristics  Is  given.  Creee  strain  Is  treated  as  an  equivalent 
external  loading.  Two  nonlinear  differential  equations  similar  to  Retainer's  equations  are 
oresentod  and  solved  by  the  iterative  scheme  and  the  finite  difference  approximation.  The 
|  multlaxiai  stress-strain-time  relationship  It  formulated  from  the  empirical  uniaxial  curve  based 
i  cn  the  second  deviated  stress  and  strain  Invariants.  Examples  are  given  for  a  corrugated  tube 
subject  to  an  axial  load  at  an  elevated  temperature.  The  results  snow  the  nonlinear  effects  duo 
to  finite  strain  and  creeo  strain.  The  method  shosvn  may  also  be  applied  to  axlsymmetrlc  shells 
sub[ecr  to  thermal  and  plastic  stralns.(  A  face-centered-njble  polycrystal  under  cyclic  loading 
:i  is  considered.  It  Is  shown  that  silo  in  the  primary  slip  system  ccuses  the  resolved  shear  stress  in 
a  second  slip  system  to  Increase  to  the  critical  value  and  then  this  second  slip  system  slides. 

Slio  In  thU  second  kilo  system  greatly  increases  the  rate  of  the  local  plastic  strain  bulld-uo 
j  in  the  primary  silo  system.  Hence  the  occurence  of  the  second  active  slip  system  may  greatly 
1  Increase  rate  of  extrusion  and  intrusion  commonly  observed  in  face-centered-cubic  metals. 
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FOREWORD 


The  research  wc.k  describee!  in  mis  report  was  performed  by  ARA,  Inc., 
West  Covina,  California,  for  the  Mechanic!  Division,  Directorate  of  Engineer¬ 
ing  Services,  Air  Force  Office  of  Scientific  Research,  Arlington,  Virginia  22209, 
under  Contract  Number  F44620-71 -C-0043.  The  research  represents  the  eleventh 
year  program  and  is  part  of  a  continuing  effort  in  the  study  of  Inelastic  theory  of 
structures  with  large  deflections  and  fatigue  damage  In  metals  due  to  cyclic 
loadings.  The  Project  Engineer  was  Dr,  Jacob  Pomerantz. 

The  studies  presented  began  1  February  1972  and  were  concluded  1  Febru¬ 
ary  1973.  Mr.  Bernard  Mazelsky  was  the  Principal  Investigator.  The  main  con¬ 
tributor  to  the  study  was  Dr.  S.  R.  Lin  of  ARA,  Inc.  Mr.  S.  C.  Liu  of  ARA,  Inc. 
participated  In  this  program.  Professor  T.  H.  Lin  of  UCLA  served  as  consultant 
on  the  analytical  study. 
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Part  A 


Large  Deflection  of  Axliymmutric  Shell  with  Creep  Strain 

I  SUMMARY 

A  numerical  method  for  analyzing  large  deflection!  of  shells  of  revolution 
under  axisymmetric  loading  with  otbltrary  creep  chcracterlitlci  It  given.  Creep 
strain  Is  treated  at  on  equivalent  external  loading.  Two  nonlinear  differential 
equations  similar  to  Relssner's  equations  are  presented  and  solved  by  the  Iterative 
scheme  and  the  finite  difference  approximation.  The  multiaxial  stress-strain-time 
relationship  Is  formulated  from  the  empirical  uniaxial  curve  based  on  the  second 
deviated  stress  and  strain  invariants.  Examples  are  given  for  a  corrugated  tube 
subject  to  an  axial  load  at  an  elevated  temperature.  The  results  show  the  non¬ 
linear  effects  due  to  finite  strain  and  creep  strain.  The  method  shown  may  also  be 
applied  to  axtsymmetric  shells  subject  to  thermal  and  plastic  strains. 

II  INTRODUCTION 

Recent  Inteiest  in  the  design  and  fabrication  of  structures  at  elevated 
temperatures  hcs  led  to  the  Investigation  of  creep  effects  in  a  number  of  structural 
configurations.  Analyses  of  beams  columns  and  plates  W  with 

creep  have  been  shown.  The  large  deflection  of  rectangular  plates  with  nonlinear 

191 

strain-hardening  creep,  was  given  by  Ho  and  Lin  1  J.  Recently  the  study  of  the 
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inelastic  oenavior  of  shells  has  become  or,  area  of  Ir, creasing  interest.  The  large 
inelastic  deformation  of  spherical  shells  has  been  studied  by  Plan  and  his 
associates  ^  The  detailed  method  of  calculating  the  stresses  and  deflect¬ 


ions  of  a  corrugated  tube  with  finite  strain  and  nonlinear  creep  hat  not  been  found 
in  the  literature.  The  present  study  develops  such  a  method  to  calculate  stresses 
at  different  points  of  an  axisymmetrlc  shell  under  axisymmetrlc  loading  such  at 


corrugated  tubes  with  arbitrary  creep  characteristics. 

Inelastic  strain  gradient  hat  been  shown  to  have  an  equivalent  effect  as 
that  caused  by  a  distributed  body  force  on  the  structure  ^  This  method  of 
equivalent  load  has  been  used  in  many  analytes  of  inelastic  bending  of 
plates  In  the  present  study,  the  same  concept  Is  applied  to  the 

analysis  of  a  corrugated  tube  undergoing  very  large  deflection  with  creep  strain. 
The  nonlinear  governing  equations  are  solved  by  an  iterative  procedure.  The 
numerical  results  demonstrate  the  neoetsary  convergence  properties  of  this  method. 


Ill  BASIC  EQUATIONS 

Using  cylindrical  coordinates  r,  B  ,  x,  the  middle  surface  of  a  shell 
of  revolution  may  be  represented  by  the  parametric  equatlonsi 


f  r  -  r<i> 

\  z  - 


O 


so  that  the  parameter  ^  together  with  the  polar  angle  &  are  coordinates 
on  the  middle  surface.  The  equation  of  the  deformed  middle  surface  Is  taken  In 
the  form 
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r  r  -  r.  i-  u 
1  r  -  z,^-  w 


(  2  ) 


where  subscript  "o"  denotes  the  quantity  referring  to  the  undeformed  state.  The 
displacements  u  and  w  are  components  In  the  'adia!  direction  r  and  axial 
direction  z  of  the  shell/  respectively, and  p  is  the  change  of  slope 
angle  9  of  the  meridians  with  respecr  to  the  z  -  axis  as  shown  In  Figure  1 . 
from  geometry/  wa  have 


cos4> 

sin  <f> 


where  S  is  the  distance  measured  along  the  meridion.  The  principal  strains  at 
the  middle  surface/  In  the  tangential  and  circumferential  directions/  are  given 
as 


(4) 


The  immedicto  consequence  of  the  definitions  of  the  components  and 
2e^  is  the  compatibility  equation 


os<f»  -<(‘•€9^')  “  -  Cnf)  (5) 
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wnere  prime  denoi.?s  ditfu.'unriaticn  with  respect  ro  ^ 
is  given  in  the  form 


and  where  & 


Or' 


o' 


Let  h,  be  the  distance  along  the  normal  to  the  deformed  shell  surface 
from  the  middle  surface,  by  using  Kirchoff’s  assumption  that  normals  to  the  middle 
surface  remain  norma!  to  the  deformed  middle  surface/  the  strain-displacement  re¬ 
lation  mcy  be  written  as 

/ 
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r  «♦- 


r  (J? - Si. ) 

>  51  c*  Ox,  ' 


l  *w,  - 


by  Hooke's  law,  we  have  the  stresses  along  principal  directions  <f>  and  0 


*“  yjp »  ( ^  -*■  J  <%/  —  ) 


<r6  -  7~rC^->  •%-<-*><; 

where  c  is  Young's  modulus/  is  Poisson's  ratio/  and/  g!  and  &* 


(i 


are  plastic  components  (Including  creep)  of  principal  strains.  From  Equations  (7) 
and  (8)/  the  stress  resultants  and  the  sectional  moments  are  found  as 
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v/nerc  r  \  ^  ».*.v  irii  of  the  shell,  or.d  the  i ."elastic  terms  with  subscript  "  I  " 

cru  cefined  as  tollows 


\  (e*+*e!>  d$ 


Consider  o  differential  element  cut  from  a  shell  of  revolution  by  two  adjacent 
meridian  planes  and  two  sections  perpendicular  to  the  meridian  as  shown  in  Figure  1 . 
l  or  the  stress  resultants  in  radial  and  axial  directions  be  denoted  by  H  and  V,  re¬ 
spectively.  The  equilibrium  conditions  give 

(rV)'-t-  cvr-ft,  »  0 

CrH  /-  <VNlfl  *-  a  r  pH  —  0  (ii) 

(r.v^ o<Mtcos<p  +  otr  ( -  I Uoi  »  o 

where  f  and  are  the  components  of  surface  load  intensity  in  the 

axial  and  radial  directions.  The  first  two  of  Equation  (11)  are  the  conditions  of 
force  equilibrium  in  tne  axial  and  radial  directions,  while  the  third  equation  is 
the  condition  of  moment  equilibrium. 

Analogous  to  corresponding  Reissner's  work  in  shells  with  finite  deflection 
but  restricted  to  small  strain  the  system  of  equations  (4),  (5),  (9)  and  (11)  can 
be  reduced  to  two  simultaneous  equations  for  the  unknown  meridian  angle  change  & 
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and  c  stress  function  ^  defined  in  terms  of  the  horizonta'  resultant  H  by 


*  =  — 
(■  Et 


'12  ) 


i he  resultant  two  differential  equations  are  of  the  following  form 
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Et 

ilex (r V) ‘ 

Ft 


-  f  [  f  aAf  t  ,>  (Siu,^  -  f W-  +  &&&  (3 
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*  •  fi  r  ' 
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«i  iiiOuiCi  go  no* co  i i *g *  1  <  1  o  1  TtciCitii  1  c  1  or m l»  ^  1  g  i  q  » i  *  c  g  •  •  •  iiio  iic^iji"*nGi*o  j< oc  o  j  *  g  g 
equation,  art;  equivalent  to  an  additional  set  or  external  loading,  in  rhis  analysis, 
no  buckling  is  considered,  hence  the  presenr  analysis  is  valid  only  before  circum¬ 
ferential  buckling  occurs. 

IV  PCLYAXiAL  STRESS- STRAIN- TIME  RELATIONSHIP  OF  CREEP 

To  calculate  the  creep  strain,  the  smooth  stress-rime  curve  has  often  been 

approximated  by  a  series  of  finite  steps,  each  of  which  consists  of  a  constant 

n  / 1 

stress  period  followed  by  an  instantaneous  Increment  of  stress  1  .  For  applying 

creep  tesr  dara  under  constant  stresses  to  the  ccsos  of  varying  stresses,  the 

rl  2' 

mechanical  equation  of  state  for  uniaxial  loading  is  generalized  to  1  J 


-  F  ( tr* 


(14) 


for  the  mulriaxial  stress  state,  where  the  effective  stress  cr  and  the  effective 


<sP 


r 

strain  <£  are  given  as.  In  the  case  of  the  axisymmetric  shell. 


_l  r 
fz 


1  ■<% 


f  V*  -  c r/  -v  cr/  J 


€  -  4  [  (<-  ep9  A  (el-  e{)l-r  ce'-  e£>* ) 


-+■  Sr. 


e'  =  0 


Equation  (15)  reduces  to 


cr'  =  <  <r#'+  a-/-  %% 
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(15) 


The  principal  creep  strain  along  the  ^  direction  has  been  denoted  by 
Since  creep  strain  produces  no  diiitation, 


(16) 


(17) 


on  r.-.c.-o-.cnrc!  time  interval  A  T,  Equation  (i  4)  and  (1  7)  give 


U  <4  A-  ( 4+  ^  4  ?+  c  *£++etx3S 

Hr-  ■*  _ 

-  -y  Fc  <r,  £  ;  -^T 

P 

7 nc  incremental  creep  strain  components,  and  A^jp  ,  are  assumed  to 


’ 1  n?  .  /rf  a ? \^1 


(IS 


oe  proportional  to  the  corresponding  deviatorlc  stress  components: 


_  *eZ  (19) 

iav-crp  2crtf-a> 

n  7] 

This  is  commonly  assumed  in  Incremental  theory  of  plasticity  .  With  known 
stresses  and  creep  strains.  Equations  (18)  and  (19)  give  the  incremental  creep 
strain  components  in  an  incremental  time  interval  Af*  .  The  total  creep  strain 
at  any  time  is  the  sum  of  the  creep  strain  increments  of  all  the  time  steps.  From 
these  creep  strains,  the  Inelastic  terms  N^  ,  N^  ,  M,^  and 

as  given  in  Equation  (10)  can  be  calculated  and  substituted  into  Equation 
(13),  to  yield  a  set  of  equivalent  loading. 


V  NUMERICAL  PROCEDURE  FOR  ANALYZING  CORRUGATED  TUBE 

We  now  consider  a  shell  with  periodic  corrugations.  Assuming  that  the 
corrugations  are  symmetric  about  their  crests  or  troughs.  One  period  of  the 
corrugations  goes  from  z  =  0  to  z  =  4c,  so  that  the  wave  length  is  4c.  Disregarding 
end  effects,  the  symmetry  of  the  shell  meridians  suggests  that  we  may  consider 
only  half  of  the  period  with  the  following  boundary  condition: 

Z  «■  O  &  2C  (  20  ) 

if  we  restrict  ourself  to  the  problem  of  corrugated  tube  subjected  to  a  constant 
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c.x;a;  i oCid  i<,  Equation  (13)  reduces  to  the  following: 


£*?-*■  =  <?, 

f  +  ^|S  ■*•  a?Y  aZ?  =  ti 


(21  ) 


wnere 


/»  1 9  ^  \  /  ^  \  i 

*3  = 

A.  «=*  —  -  — * 

**  r  <x. 

*«-  "(•£ ' au+.  +  Olismt.  )  ( 

«,= 

^  “  ~r  {.*“%'* 4>X  -  $"40*4  -  •£ 

+  -p- (0tocoZ# - <£'- _  ILHzlL**. Sir><f>  ] 

€  —  -TrC -7-' */«<?»  ■+)J(’~Sih</>  -  -^rf/rxp  •*-<t>4>£<t>)] — (!— 

C It  1  *•  '•  *  21T£t 

”  ~ J{ "*  (  )Coit^  -4-  (Vn^-  p  >  Stn  <fc  J 

+T^‘  ^ ;  -  "h  *•■  ;]  -  <^  (~  -  *^  -fp) 

^  "  “r  [•^2^>4>'  -*•  <?“*»<}>  H t  -*>*>*?  -Scoift  -  ^ 


(22) 
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Equation  (21)  is  solved  by  finite  difference  and  iterative  technique  as  shown  in  the 
following. 

On  the  basis  of  conventional  central  difference  formulae 


tne  differentia!  Equation  (21 )  can  be  written  in  finite  difference  form  at  any 
stations  n  except  at  z=o  and  2c, 


(24) 
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(28) 


Tne  set  of  equations.  Equation  (26)^111  be  solved  by  the  procedure  equivalent  to 
solution  by  the  method  of  Gaussian  elimination  used  in  other  problems  of  shells 

ns  1 91 

of  revolution  '  J.  The  elimination  technique  proceeds  as  follows:  Let  rhe  two 
points  z=o,  2c  be  denoted  by  n=o  and  n=N  respectively,  so  that 


(3  =  Q  «  0 


(29) 


From  symmetry,  we  also  have 

f  =  "  "  Pm-, 

1  7-,  "  "  7.  ;  =  -  7m-. 


(30) 


fhe  first  of  Equation  (26)  i.e.  n=l  is  solved  for  y?  in  terms  ^  and  : 
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f,  “  (t'-Xfi-a:?,) 

rs 


(31  ) 


This  result  is  substituted  into  the  next  equation  n=2  and  pz  will  be  found  in 
rerms  of  Pi  and  ,  and  so  on.  Finail/  the  last  equation  of  Equation 

(26)  will  determine  7 .  ,  and  then  with  7,-0  >  all  of  the  7  s  will 

c-e  calculated  in  reverse  order.  Accordingly,  we  may  write  as 


?n  =  R„-  P»/V,  -  <2 n?, 


n+i 


(32) 


:rom  Equations  (29)  through  (32),  we  have 


4 


f  *  ~-sh> 


K  =  o 


s 

/ 


p  =.  _ gj 

BJ-  <V 


;  p.  -  0 

L  ~  —  0 

Substitution  of  Equation  (32)  into  Equation  (26)  yields  the  result: 


_  3/ 


(33) 


where 


Rn  = 


1  Pn 

l  a. 


i£  K-,  -  £/  <  R,.2  -  P„„  rJ 


p  ^  @n-i  +  Of  Pn-i  Gln-t  7 


IZ. 

*>n 


(34) 


n=*  z,  •••,  v-2 


=  Bjn-  BP  P*_,  -r  Bsn  CP„.,  -  Qn„) 


i  m 


Is  recurrance  relation  together  with  the  initial  values  from  Equation  (33)  give 
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o!i  rhe 


R  's,  P  's  ana  Q  ‘s  up  ro  R  ,  P  and  Q  .  I  he  last  of  Equation  (26),  r.-N- 
n  '  n  n  r  n-2  n-2  n-2 


and  symmetry  condition.  Equation  (30),  determine  ^-\  as 

ft  _ 


(35) 


wnere  Rfi_,  and  ^  are  the  extension  of  Equation  (34)  to  n=N-l .  With  ^ 

ond  ?>i->  known'  h- i  '  /V3  ,  "  •  .?*.>?,  can  be  found 

from  Equation  (32)  and  then  ell  the  ^  s  can  recdiiy  be  obtained  from  Equation 
(27).  This  process  is  readily  programmed  for  the  high  speed  digital  computer. 

The  stress  resultants,  meridian  strains  and  displacements  defined  in  accordance 
wirh  Figure  1  are  obtained  from  p  and  ^  ■  as  follows 

(  y  COS  -i-  sin  <P  ) 


—  -fr  (- ?  Sir>4>  -t-  ~^:CoiCp  ) 

~  f  f'V  *♦*  ^  (  1^3  ■+•  i^Q£  )  ] 

=  ~Et  ^  ,Nls  *  ^  ] 


(36) 


r° 


—  f  *.[</*  ^^0  ^  —  Sm  <$0  J 


In  the  Inelastic  problem,  the  loading  path  Is  important.  The  time  is 
increased  in  Increments  In  the  case  of  creep  deformation.  In  each  step  of  time 
increment,  the  initial  values  of  T  f  and  Bq  at  each 
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station  arc  known.  The  stresses  ,  U"^  ,  CT 


and  the  strains  ,  €. 


** 


ccr.  oo  caicu.ctea  from  equations  (4),  it  (3y  ana  \17).  *ne  incremental  creep 


strains  A  and  Ad?/  in  an  incremental  time  interval  aT  are 
obtained  from  Equations  (13)  ana  (19)  with  given  Function  F  (  CT  t  &  )• 

inelastic  stress  resultants  and  sectional  moments  are  found  from  Equation  (101.  With 
assumed  values  of  6  and  p  ,  the  coefficients  A;  s  and  A;  s 
can  be  calculated.  Equation  (24)  can  then  be  solved  for  new  values  of  ^ 
and  p  .  The  procedure  is  repeated  until  the  differences  between  the 
successive  p>(s  and  ps  are  within  desired  tolerances.  The  calculation  can 
then  be  extended  to  the  next  time  increment,  in  the  following  numerical  example, 
the  tolerance  for  of  0.01%,  was  used. 

Vi  NUMERICAL  EXAMPLES 

The  above  procedure  is  first  applied  to  the  static  deformation  of  a  corrugated 
tube  similar  to  the  one  tested  by  Donnell  ^  .  The  meridian  consists  of  circular  arcs. 
The  dimensions  (Figure  2)  are 
f  a  =  5 . 33  I  n 
b  =  0.235  in 
c  =  0.670  in 
^  t  -  0.065  in 

and  E  =  30x10^  psi,  O  =0.3.  Figure  3  shows  the  axial  load  deflection  curve 

for  tensile  and  compressive  loadings.  The  linear  solution  checks  very  well  with 

1211 

rnat  of  Ciark  and  Reissner  1  .  It  is  seen  that  the  corrugated  tube  under  tensile 

axlai  force  is  stiffer  than  under  compressive  axial  force.  This  was  pointed  out  by 
*22 1 

homada  1  ‘ t  in  the  elastic  large  deflection  theory  of  shell  of  revolution. 

The  incremental  procedure  is  also  applied  to  a  7075-T6  aluminum  alloy 
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coivu£c;i«2ci  'obis  undergoing  large  doformcrion  ar  600°  F.  Tlio  geometric  dimensions 
are  same  as  above.  The  cpplied  axial  load  is  2124  ibs.  The  uniaxial  creep 
characteristics  of  the  material  are  approximated  by  the  following  stress-strain-time- 
ru!  a  t;  on: 

*.  **  JL-  +  A  (3<r) 

E 

where 

A  5.25  x  10"7/hr 
B‘»1.92x10*3in2/lb 
k  ■  0. 66 
E  =  5.2x  106  pel 
e  *  total  (train .. 

(T  ■  stress  In  psi 
T  “  time  In  hour 

and  Poisson's  ratio  j/  is  0.32.  The  creep  characteristics  are  assumed  to  be  the 
same  In  tension  and  In  compression.  The  first  time  Increment  was  taken  to  be 
0.0001  hr.  The  subsequent  time  increment  was  increased  In  each  step  at  constant 
rate.  There  were  100  time  Increments  to  reach  the  total  of  16.66  hours.  The 
100th  time  Increment  was  .5050  hr.  The  number  of  finite  difference  stations  along 
the  midsurface  in  half  period  of  corrugation  was  49.  The  number  of  grid  points 
ucioss  the  thickness  was  13.  The  tolerance  used  In  the  presor.t  calculation  was 
0.01%  for  *1  .  To  obtain  this  degree  of  accuracy,  3  cycles  of  Iteration  was 

generally  found  to  be  sufficient.  The  computation  was  carried  out  onan  IBM  360 
and  the  computing  time  for  this  creep  analysis  was  within  one  minute,  figure  4 
shows  the  Increase  of  total  axial  deflection  with  time.  The  variations  of  the  extreme 
fiber  stresses  at  two  ends  at  different  instants  are  shown  In  Figures  5  and  6. 
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v;:  conclusions 

A  numerical  method  is  shown  for  the  creep  analysis  of  axisyrnmetric  shells 
subject  to  cxisymmetric  loadings  with  large  deflections.  A  numerical  example 
of  creep  deformation  of  a  corrugated  tube  under  an  axial  load  with  finite 
deflection  is  given.  Vertical  aisplacement  and  stress  distributions  of  the  tube 
at  different  instants  of  time  cfter  loading  are  shown.  This  method  may  be 
extended  to  analyzing  elastic-plastic  deformation  of  shells  of  revolution  with 
finite  deflection. 

In  the  numerical  example  of  creep  analysis,  ICO  time  increments  were  used. 
The  time  increment  increases  from  an  initial  value  of  0.0001  hr.  to  0.5050  hr.  for 
the  100th  increment.  The  convergence  of  the  solution  with  the  number  of  finite 
difference  stations  along  the  meridian  was  found  to  be  good.  It  was  found  that  the 
use  of  25  and  49  finite  difference  stations  gave  only  0.2%  difference  in  total  axial 
deflection.  The  number  of  grids  across  the  thickness  and  time  increments  used  also 
provided  good  convergence. 
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Part  B 


Effect  of  the  Secondary  Slip  System  on  Early  Fatigue  Damage 


I  SUMMARY 

A  face-centered-cubic  polycrystal  under  cyclic  loading  is  considered. 

It  is  shown  that  slip  in  the  primary  slip  system  causes  the  resolved  shear  stress  in  a 
second  slip  system  to  increase  to  the  critical  value  and  then  this  second  slip  system 
Siiaes.  Slip  in  this  second  slip  system  greatly  increases  the  rate  of  the  local  plastic 
strain  build-up  in  the  primary  slip  system.  Hence  the  occurrence  of  the  second 
active  slip  system  may  greatly  increase  rate  of  extrusion  and  intrusion  commonly 
observed  in  face-centered-cubic  metals. 
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1  INTRODUCTION 

Fatigue  failure  is  characterized  by  fracture  under  repeated  stress  far 
below  the  static  fracture  strength  and  usually  occurs  by  the  processes  of  crack 
nucleation,  propagation  and  final  rupture.  Even  though  regions  near  the  nucleus 
of  the  crack  and  along  its  path  of  propagation  may  be  highly  distorted,  material 
away  from  these  regions  may  not  show  any  significant  macroscopic  plastic  deforma¬ 
tion.  In  a  single  crystal,  the  fatigue  crack  often  lies  along  a  slip  plane.  Ewing 
and  Humphrey  ,  studied  the  formation  of  cracks  by  the  microscopic  observations 
of  the  surface  appearance  of  specimens.  Their  tests  on  Swedish  iron  revealed  that 
very  few  slip  lines  were  produced  by  the  first  few  thousand  cycles  but  as  the  test 
progressed,  more  slip  lines  formed.  These  slip  lines  are  the  intersections  of  sliding 
planes  with  the  free  surface.  As  stressing  continues,  bands  of  slip  are  produced 
and  intensified.  Slip  lines  in  fatigue  specimens  were  found  to  be  spaced  not  as 
evenly  as  those  in  unidirectional  loading.  Fatigue  cracks  eventually  are  formed 
in  these  intense  slip  bands. 

The  first  quantitative  physical  theory  of  fatigue  crack  nucleation  based 
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on  the  microstress  field  caused  by  slip  was  proposed  by  Lin  and  Ito  1  ' 

The  theory  was  applied  to  study  the  influence  of  strain- hardening,  grain  size  and 
mean  stress  on  early  fatigue  damage  1  '  \  The  results  show 

that  the  increase  of  the  strain- hardening  rate  and/or  the  decrease  of  the  grain 
size  decreases  the  rate  of  early  fatigue  damage,  and  the  effect  of  mean  stress  on 
the  number  of  cycles  to  failure  was  found  to  lie  between  the  values  predicted  by 
Gerber's  parabola  and  by  the  modified  Goodman's  line. 

The  Lin-lto  theory  was  in  part  based  on  the  experimentally  observed 
formation  of  slip  bands  and  is  briefly  reviewed  here  for  reference.  The  slip  lines 
formed  in  forward  loading  and  those  formed  in  reverse  loading  are  very  closely 
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spaced  but  dlsvincr  from  each  other  and  occur  Sr.  the  same  slip  band.  Positive  slip 


occurs  in  one  siip  line  In  forward  loading  and  negative  slip  occurs  in  a  nearby  slip 
line  In  reverse  loading.  Since  lattice  imperfections  exist  In  all  metals.  Initial 
heterogeneous  microstress  fields  exist  Inal!  metals.  The  initial  stress  field  favor¬ 
able  ro  this  type  of  siip  formation  is  clearly  one  producing  resolved  shear  stresses 
of  opposite  sign  in  two  closely  spaced,  parallel  thin  slices.  As  shown  in  Figure  7, 
two  thin  slices  in  a  most-favorably-oriented  crystal  located  at  the  free  surface 
have  an  initiol  resolved  shear  stress  field  which  Is  positive  in  slice  P  and  negative 
in  slice  Q  .  Under  tensile  loading,  the  positive  critical  shear  stress  Is  attained 
and  hence  siip  occurs  in  P  .  Since  the  residual  stress  field  caused  by  slip  is  con¬ 
tinuous,  slip  in  P  relieves  the  positive  resolved  sheer  stress  not  only  in  P  but 
also  in  Q  .  This  helps  to  keep  Q  from  reaching  the  positive  critical  shear  stress 
in  tensile  loading,  and  also  increases  the  negative  resolved  shear  stress  that  causes 
Q  to  slide  more  easily  in  reverse  loading.  When  the  negative  resolved  shear 
stress  in  Q  reaches  the  critical  shear  stress  in  reverse  loading,  slip  occurs  in  Q 
and  a  r.ew  residual  stress  field  is  produced.  The  new  stress  field  increa^s  the 
positive  resoivea  shear  stress  not  only  In  Q  but  also  in  P  .  This  increase  of 
positive  resolved  shear  stress  in  P  causes  P  to  slide  more  easily  in  the  next 
tensile  loading.  This  process  is  repeated  for  every  cycle  and  builds  up  the  local 
plastic  strain  at  these  two  thin  slices,  P  and  Q  .  The  plastic  strain  in  P  and 
G  vends  to  push  out  the  region  between  them  and  starts  an  extrusion.  If  the  signs 
o;  Ir.irial  stress  in  P  and  Q  are  changed,  an  intrusion  instead  of  extrusion  wiil 
jc  •  rii  1 1  ated.  Thus,  the  continuity  of  the  relief  of  stress  provides  a  natural  gating 
mecr.a.'"  m  for  the  monotonic  build-up  of  local  plastic  strain  and  no  special  gating 
mechanism,  such  as  a  specific  movement  of  the  Lomer- Cottrell  barrier,  as  proposed 

T23 1 

u/  Kennedy  1  J  is  necessary. 
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in  Lin  and  iro's  ca.cu.ation,  ir  was  ossumec  the.;  oniy  the  mosr-favorably- 


orien.ed  slip  system  would  be  activated.  This  system  is  referred  to  as  the  primary 
s.ip  system.  For  face-cer.rered-cubic  (fee)  metals  such  as  aluminum  and  its  alloys, 
each  crystal  has  four  slip  planes  with  three  slip  airections  on  each  plane  and  hence 
there  are  twelve  slip  systems.  During  the  process  of  extrusion,  the  material  is 
elongated  along  the  direction  of  the  extrusion  and  hence  a  tensile  stress  is  developed 
in  the  region  between  ?  and  Q.  Similarly,  a  compressive  stress  is  developed 
ciong  the  direction  of  an  intrusion.  This  stress  will  cause  other  slip  systems  to 
slide.  Tne  purpose  of  the  present  study  is  to  evaluate  the  slip  in  these  secondary 
slip  systems  during  the  early  stages  of  fatigue, 
iii  RESOLVED  SHEAR  STRESS  FIELD  CAUSED  BY  SLIP 


Among  the  different  mechanisms  of  plastic  deformation  in  metals, 

crystallographic  slip  has  been  shown  to  be  the  principal  process  of  plastic  deforma- 

129 1 

rion  ir.  fcce-centered-cubic  metals  at  low  and  intermediate  temperatures  *.  From 
single-crystal  tests,  it  has  been  shown  1  '  J  that  slip  depends  on  the  resolved  shear 

stress  and  is  independent  of  the  normal  stress  on  the  slip  plane.  The  resolved  shear 
stress  necesscry  to  initiate  or  to  cause  further  slip  is  generally  known  as  the  critical 
shear  stress.  Those  slip  systems  in  the  crystal  with  highest  resolved  shear  stress 
are  called  the  most-favorabiy-oriented  slip  systems  and  are  the  first  to  reach  the 
critical  shear  stress.  When  the  critical  shear  stress  is  reached  in  some  region  in  the 
most-favorabiy-oriented  crystal  in  the  body,  slip  occurs  and  causes  plastic  strain. 
Should  the  load  be  removed  after  sliding  in  the  region,  the  slip  would  remain  and 
cause  a  residual  stress  field.  If  the  load  is  applied  again,  assuming  no  further  slip 
during  this  reloading,  the  residual  stress  plus  the  applied  aggregate  stress,  gives 
the  resultant  stress  field  in  the  body.  To  find  the  residual  stress  due  to  localized 
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s.ip,  fhe  analogy  between  the  plastic  strain  gradient  and  a  body  force  In  causing 

stresses  in  an  elastic  medium  first  shown  by  Reissner  and  developed  by 

Esheiby^^  andLin^1^  Is  used.  Referring  to  a  set  of  rectangular  axes,  let 

the  plastic  strain  tensor  due  to  slip  be  denoted  by  €■■  .  The  strain  distribution 

£ 

in  a  body  with  plastic  strain  under  external  load  is  the  same  as  that  in  a  purely 
elastic  body  with  the  additional  equivalent  body  force  and  surface  fraction: 

f  F;  =  -  tpe-.y, 

where  JA.  is  fhe  shear  modulus,  P  is  the  component  of  the  unit  surface  normal  P 
on  x.  axis,  and  subscript  after  a  comma  denotes  differentiation  with  respect  to 
a  coordinate  variable.  The  residual  stress  field  of  the  actual  body  will  be  given 
as: 


0) 


% 


Zn  -  ZVL<2: 


(2) 


where  Z  -  is  the  stress  field  caused  by  F-  and  T.  . 

Js  1  i  i 

The  aggregate  considered  is  of  fine  grain;  grain  size  being  small  com¬ 
pared  with  fhe  total  volume  of  the  aggregate.  On  Figure  7,  x^ ,  end  Xg  are  a 
set  of  orthogonal  cartesian  coordinates,  where  the  axis  is  parallel  to  the  width 
of  slices  P  and  Q  ,  and  x^  axis  is  normal  to  the  free  surface.  Since  the  two 
slices  P  and  Q  are  parallel  and  the  thickness  of  each  slice  is  much  less  than 
its  width,  the  deformation  state  can  be  considered  under  plane  strain.  By  using 
Airy's  stress  function,  the  stress  distribution  caused  by  a  point  force  in  a  semi- 


infinite  plate  for  plane  stress  has  been  givan  by  Melon 


[34] 


Melon's  solution 


[35] 

was  modified  for  plane  strain  by  Tung  and  Lin  .  This  stress  function  is  used 
in  rhe  calculation  of  slip  strains  and  stresses  in  polycrystalline  aggregates  near  the 
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vree  surface.  Let  Z-.  be  the  stress  at  X  “  (  X, ,  X2  >  due  to  a  unit 

'i 

force  acting  at  X  a  (  X,  ,  along  the  direction  of  x^  axis,  and 

^(x,*;  be  the  corresponding  stress  function;  we  have: 


Z  lx  x)  =  -2^- 

where  the  stress  functions  are  given  as 


K  «-  1,2 


x,x,(x,t-yfj 


4><X,X>  —■  «*•  ^p— 

[  $j0c,O-  +  -  2p^~ 


r  A.TU-VJ 
o0  =s  pa-zc-; 


$t  *=*  "han 


-I  x.-x. 


-fan 


*»  Xj-Xj 


-  IT  «  S,  <  *C 
T  <  I 


|  Z,-  (XrX,^V  <X,-Xa> 

^  Z*=  <x(-4-x,/f  (Xw~x»/ 

It  should  be  pointed  out  that  Equations  (3)  and  (4)  are  valid  only  when  XAf  X  . 

Consider  a  point  x  in  a  crystal  that  has  slid.  Let  m)  be  normal 

to  the  slip  plane  and  J  Cm)  be  In  the  sliding  direction  of  the  m'  active  slip 
system  at  the  source  point  x,  and  ^(x,  m)  be  the  plastic  shear  strain  caused 
by  this  slip.  This  plastic  strain,  referring  to  x -coordinates,  is  then 


e“  (  x  >  <=»  rn.. 

li  't 
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Wi»vm*£ 


m 


ij.  -  <rrO  £  (m  > -<-  ^,<rn; 


In  which  £.(m)  end  (nr.)  denote  the  components  of  unit  vectors  along 
;•  (m)  and  (m)  directions,  respective! y.  It  is  shown  in  the  appendix  that 
I  r..<x,jo  *— = — — civ,  </£,  converges  uniformly  in  the  slid  region  S  ,  hence  by 

J j.  'i  &x| 

equations  (1)  and  (6),  the  stress  at  field  point  x  caused  by  a  given  plastic 
strain  distribution  e"(x, m)  in  the  slid  region  S  is  given  as 


(7) 


v-  - 


-  2 


r.V  JO  'left  dXi 

•J  2>  X* 


(8) 


This  Integra!,  and  hence  T(--  in  the  above  equation  exists  and  represents  a 

*s 

continuous  function  on  S  .  In  equation  (8),  the  repetition  of  slip  system  m 
denotes  summation  over  ali  sliding  slip  systems  at  point  x  while  the  subscripts 
k  and  H  range  from  1  to  2.  The  above  discussion  holds  also  when  the  slid 
region  S  Is  a  multi- connected  region  where  the  slip  strain  is  continuous  inside 
each  simply-connected  subregion. 


Slip  in  any  region  depends  on  its  resolved  shear  stress.  The  resol 


vea 


th 

shear  stress  in  the  n  slip  system  at  a  field  point  x,  caused  by  slip  throughout 

rhe  slid  region  S,  is  obtained  by  transforming  the  stress  "ZIg  <jo  given  by 

equation  (8)  from  the  (x^ ,  X2)  coordinates  to  the  ^  )n  coordinates,  where 

•j  (n)  Is  normal  to  the  slip  plane  and  ^  (n)  Is  along  the  slip  direction  of  the 
th  . . 

n  sup  system  at  x: 


~Ce  iXj  r\> 


2  ( 


(9) 


in  wnlcn  n. .  is  defined  by  Equation  (7),  The  substitution  of  Equation  (8)  into 


ll 
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Equation  (9)  yielas 

f  K 

Z^(X,o.)=  -p.  rn^  n.j  j  Zjj  (x^xi  dxtdx1 

s 

-  ,  K 

Letxlng  (x,  n;  x,  rr«)  =  ~  n. .  mj<jj  2Tjj  £x,  Jo  /  *he  above  equation  may  be 
written  In  the  following  form. 


—  -z/tj  a J}(Ko.  x. 


*xi 


It  is  shown  in  the  appendix  that  the  integral  in  Equation  (A-3)  over  the 
region  Sf(x)  tends  to  zero  as  the  radius  £  of  S£cx)  approaches  zero.  Hence, 
Equation  (1 1)  can  be  written  as: 


Z  <x,n>  = 


-iu  f  A L(x,n.  x,rr>'j  Vx,  Jx^  +  0  ) 

1  s-s,  J  **< 


as  £ 


Since  the  integrand  in  the  above  equation  is  regular  in  the  integration  domain 
S  -  S.  tne  area  integral  can  be  integrated  by  parts  and  particily  transformed  into 
a  line  Integra!  by  means  of  Green's  lemma: 

Z^cx^n}  **-zu(p  {  X,(xtn,  x,m>e"<x,r/i)tlxl-  hl2(x,n.  x,/n>e'(x,m)dx,  J 

+  ^  *.m)  -*■  eCx,^)  dx,dxx  03] 

S~Sg 

-*■  OCcr>^  as  £  O 

whore  P  and  17  ,  taking  counterclockwise  as  positive,  denote  surfaces  bounding  the 

regions  S  and  Sg  respectively.  The  slip  strain  in  each  region  of  a  given  slip 

orientation  is  taken  to  be  continuous  inside  and  drops  to  zero  at  the  boundary.  The 

line  integral  In  Equation  (13)  along  P  vanishes.  For  numerical  calculation,  the 

slid  region  is  divided  into  thin  parallelogram  regions  unless  it  is  otherwise  specified. 

To  simplify  calculations,  plastic  strain  in  each  grid  is  assumed  to  be  constant.  Let 

S,  denote' the  region  of  f  1  grid,  fj.  the  plane  surface  bounding  the  grid  and 
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jo  the  grid  containing  Sc  .  Then  Equation  (13)  becomes 

'csi*‘n>  ~  Zj*£  (Xt.mlS  [  Ni ,<*,n,  x,rr<)dx>.-  sijOqn,  x,m)clxt  j 

x,my]dx,dxt 


0(Z^' . 


as 


th 


in  the  above ,  e"(x.  m)  represents  the  plastic  strain  due  to  slip  in  the  m  slip 
1  / 

system  of  the  grid  ,  wheie 

ST  *=•  S,  “  Se  H<x,  *o  > 


(14) 


(15) 


In  which  the  two-dimensionai  step  unit  function  H  (x,  xq)  Is  defined  as 


f  H  (x,  xQ)  ~  | 


if  x  £  S9 
if  x  %  S0  . 

fhe  repetition  of  r  and  m  denotes  summation  over  all  grids  S  and  all  sliding 


H  (x,  xq)  =  q 


slip  systems  at  grid  centroid  xf.  Since  the  area  Integrals  over  the  regions  Sr  in 
Equation  (14)  are  proper,  Green's  lemma  can  be  applied  to  reduce  area  integrals 
to  line  integrals.  The  resulting  equation  becomes 

n * 

it  is  clearly  seen  that  the  above  Integral  represents  the  resolved  shear  stress  at  x 
i  th 

in  the  n‘  slip  system  due  to  constant  unit  plastic  strain  caused  by  slip  in  the 

^  |» 

m  ‘  slip  system  in  the  region  S^., 

If  more  than  one  slip  system  slides  at  field  point  x,  the  plastic 
resolved  shear  strain  at  x  in  fhe  n  slip  system  due  to  slip  in  different  slip 
systems  at  fho  same  point  x  is  given  as: 


(16) 


£*(X,  -  4  nijrr)lj  & 


07) 
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'.7.o  rc7wuQi  resolved  shear  stress  caused  by  plastic  strain  in  the  site  region  S  is 


rR<x,n>=*  Zs<x,n>  -  2fLe*(x,n) 

The  substitution  of  Equations  (16)  and  (17)  Into  the  above  equation  yields: 

ZTsoqn.>  =  -  zj<  e.\xrtm  )  <5(x,n.  xr,m) 


wnere 

6utn.xr,m)  “  n-m..  H', xftm)dxt) 

'  JVr 

The  function  G(x,  n.x^,  rr.)  is  called  the  stress  influence  coefficient  of  the 

residual  resolved  shear  stress  at  x  in  the  r?n  slip  system  due  to  a  constant  unit 

plastic  strain  distribution  In  the  subparalleiogram  S,  with  centroid  at  x  ,  caused 
th 

by  o. ip  in  the  m‘  slip  system. 

The  poiycrystal  is  subjected  to  an  alternate  tension  and  compression 
.oabir.g  along  the  X2“axis  as  shown  in  Figure?.  Since  the  maximum  shear 
srress  occurs  in  a  plane  at  4 5°  with  the  direction  of  loading,  the  primary  slip 
direction  and  the  normal  to  the  primary  slip  plane,  referred  to  as  the  first  slip 
system  of  the  most  favorably  oriented  crystal  at  the  free  surface  are  at  45°  with 
the  X2  -axis: 

f  ?(ii  “  <-j? .  5?.  °> 

*•  )“>  “  <  .  Tr  ,  0  ) 

it  has  been  pointed  our  that  under  cyclic  loading,  two  distinct  closely-spaced 
slices  slide.  One  slice  slides  during  the  forward  loading,  the  neighboring  slice 
slices  during  the  reversed  loading.  The  build-up  of  large  local  plastic  shear  strain 
in  the  primary  slip  system  in  these  two  thin  slices,  as  shown  by  Lin  and  Ito  ^4, 25]^ 
tunas  to  start  an  extrusion  or  an  intrusion  in  tho  region  between  these  slices.  Con- 
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sequonrly,  an  appreciable  direct  stress/  tension  tor  extrusion  and  compression  for 
Inrrusion,  will  occur  In  this  region  between  these  slices  along  the  primary  slip 
oirecrion  (i ).  This  direct  stress  gives  maximum  shear  stress  along  x  (  direction 
on  the  x 2  piane,  or  along  ^  direction  on  the  Xj  plane.  A  f.  c.c.  crystal  has 
twelve  slip  systems.  This  direct  stress  causes  different  resolved  shear  stresses  in 
rhe  different  slip  systems.  The  one  with  the  highest  resolved  shear  stress  caused  by 
this  direct  stress  is  referred  to  as  the  secondary  slip  system/  which  should  have  the 
ill p  direction  and  the  normal  to  the  slip  plane  close  to  x-f  and  ^  axes.  To  simplify 
the  numerical  calculation/  this  secondary  slip  system  is  assumed  to  have  the  slip 
direction  and  the  normal  to  the  slip  plane  along  x-  ana  ^  axes,  i.e., 

f  ?<!J  -  < 

C  -  <:  i  ,o,o) 

The  grids  are  taken  to  be  rectangles  for  the  secondary  slip  system  in  the  x^-X2 
plane,  as  shown  in  Figure  8. 

To  find  the  stress  influence  coefficient  G(x,  n;  xf,  m),  the  surface 
integral  in  equation  (16)  has  to  be  evaluated.  The  detail  integration  of  the  inte- 

i  __  r3£j 

gral  Is  given  by  Lin  1  .  The  resulting  expressions  are  given  as  follows: 


r  4<x./;  Xr.tj  «  -  K'CK.x,.  Vr.Vrt) 

i 

!  -  '?> 

i 

■  <J7  <X,2;  X(t  l  )  «=S  “JZiJXi.Xtj  X,r+£  t 

(  Xy-f  ,  Vi  +  'T>5  “  <*'.**>  X,~T.  V*  *7 ■> 

j  %r,  ^  “  ~  (  A.  Kr-*-  ~  ,  ■+  (Xi.Xt.,  .  *>f~  §0 

|  +  (*•.*>>  V7-  <*>.*; 

v,  ^  ^  ^  ~  ^  f  x,f-+  ~ ,  x^-  li  ^ 

t-  JrJ; 


(21) 
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wnere 


--  ct/v.  v  r  s  _  ,1  f  j _ -i  ^  u(u+&)  U+2UV- </+/</«/ .  w-n/;f«-/j  ^ 

f  -jj-t-kw  V*--^ - is; - -*- - xT~  J 


» 

l 


*,.*,>=  f  [t%-|  -  wioagb " j 

2,  *u>  -  ?  ( -f-'  -  ^4^4  j 


(22) 


in  which 


V  W  *  X,  -  X, 

V/  “  X,*X, 

(  w  a  X,-X, 

2,  =  wV  ixJ  1 

a.  & 

V  Xv  *»/■+■  u/ 


(23) 


ana,  a,  o,  #  and  p  are  the  width  and  thickness  of  the  primary  grid  and  second¬ 
ary  grid/  respectively. 

When  the  thickness  of  grids  is  much  smaller  than  its  width,  i.e.,  b<c£  a, 
it  can  be  shown  that 

4<Kij  yr,  i)  zb  [J£m( v4'  x*,+t>  -X (*,. x* .  v f ,  (24) 


wncre 


K.iX'.x*.  x3  j 


•  ..  ^  i  1  4  1 

^  ^  £t/  (u-'J)  _  ,>  ✓  2u/'hUM-j-2Wn/-\/<*J-S\/  4\/(u~</  jOZ+ufJ  i 

'  *■"’  2,  X*  *  **  ^  J 


(25) 


<f,is  dearly  shows  the  linear  dependency  of  the  relief  of  resolved  shear  stress  on 
/he  plastic  strain  and  the  grid  thickness  “  b  ",  The  same  result  has  been  pointed  out 

'37 1 

by  Lin  and  Iro  .  However,  the  presen*  expression  of  stress  influence  co- 
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efficient  is  in  g  much  simpler  form  then  they  previously  given  by  them. 

.V  LOCAL  PLASTIC  STRAIN  BUILD-UP  UNDER  CYCLIC  LOADING 

The  maximum  local  plastic  strain  build-up  at  the  free  surface  is  used  as 
the  criterion  for  estimating  early  fatigue  demage.  To  determine  quantitatively  this 
plastic  strain  build-up  for  the  fatigue  nucleation  mechanism/  the  following  proce¬ 
dure  is  used. 

Consider  a  polycrystalline  metal  with  initial  resolved  shear  stress  Z*£. 

This  metal  is  subjected  to  an  alternate  stress  2T-,  .  Before  slip  occurs/  the  applied 

9 

resolved  shear  stress  2^  is  uniform.  The  resolved  shear  stress  at  x  in  the  n 
slip  system  due  to  the  applied  cyclic  stress  is 


~  n;-  Zr  “■  L(xtni  Z-,>  (26) 

where  L  (x,  n)  is  the  Schmid  factor  which  gives  the  ratio  of  the  resolved  shear 
stress  in  a  particular  slip  system  to  the  applied  load  . 

Under  the  initial  forward  loading/  the  resolved  shear  stress  is  Z**2T-  . 

When  the  critical  shear  stress  is  reached/  plastic  strain  occurs  and  produces  residual 
resolved  shear  stress  as  given  by  Equation  (19).  The  total  resolved  shear  stress  in  the 
n‘  slip  system  at  x  is  the  sum  of  the  applied,  initial  and  residual  resolved  shear  stresses 

Z<x.nj=  Lu.nJ  £.  ZxU,nj-  LjA.<q(x,nJ  (27) 

The  repetition  of  m  and  r  denotes  summation  over  all  sliding  slip  systems  in  all 
slid  regions.  The  initial  stress  field  does  not  change  with  loading,  so  that 
•— — ~ — ~  =  0.  The  incremental  rate  of  resolved  shear  stress  with  respect  to 
applied  load  is 


cl  C(K,r))  ,  ,  s  de  <X,m) 

— : ■  —  —  Ltx.m  -  iu.6m.0j  — — 

oi  c.  /  ** 
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Those  grids  in  which  rhe  magnitude  of  the  resolved  shear  stress  is  less 
than  the  critical  shear  stress  "Cc  are  referred  to  as  non-active  slip  grids.  If  Rq 
Is  one  of  these  slip  grids, 

|  ro<p,roj<  ^  /  =  0 

where  x  is  the  centroid  of  the  grid  R  .  In  the  region  currently  sliding,  the 
P  P 

magnitude  of  the  resolved  shear  stress  equals  the  critical  shear  stress.  Let  R<-  be 
one  of  the  grids  in  the  sliding  region  with  centroid  at  xg.  It  follows 

J  £  (x,,  f)>  J  =.  rc  ,  j  j  ^  o 


where  n  denotes  the  sliding  slip  system  at  xs- 

The  macroscopic  plastic  strain  of  the  crystal  represents  the  overage 

n  21 

value  in  the  crystal  1  .  Since  the  slip  is  highly  concentrated  in  the  thin 

slices,  the  macroscopic  plastic  strain  is  much  less  than  the  local  plastic  strain  ^4,25] 
The  rate  of  s  tr  a  i  n- h  a  rd  en  i  n  g  in  terms  of  the  local  plastic 
strain  is  hence  much  less  than  that  in  terms  of  macroscopic  strain.  To  simplify  the 
calculations,  this  local  strain- hardening  is  neglected.  Then  Equation  (28)  gives 

ZJA  -j  z  = 


where  x^  and  xr  denote  the  centroids  of  sliding  regions.  This  is  a  set  of  linear 
ecuations  with  as  many  non-zero  unknowns  as  there  are 

eejarions.  The  plastic  strain  increment  4e"  (xr,m)  in  the  sliding  grids  for  an 
increment  of  applied  stress  &Ta  can  readily  be  determined  from  the  value  of 


clZ. 


Substitution  of  (\rn}  into  Equation  (28)  yields  the  rate  of 


cnange  of  the  resolved  shear  stress  at  all  points.  From  the  known  values  of  £(x,  n) 

ar  rhe  non-active  points  and  the  corresponding  values  of  ,  the  increments 

c(  Z0 


I r,  'c0  required  for  each  no n- a c five  point  to  Initiate  slip  can  be  calculated  and 
compared.  The  minimum  or  these  required  increments  in  "Za  Is  applied,  resulting 
in  one  additional  active  point  for  the  next  load  interval.  During  reversed  loading, 
the  incremental  aggregate  stress  in  reverse  loading  required  for  a  new  slip  region 
to  slide  can  similar!/  be  found. 

It  should  be  noted  that  the  incremental  calculation  of  slip  strains  and 
stresses  under  quasi-static  cyclic  loading  outlined  in  the  foregoing  is  rigorous  with¬ 
in  the  framework  of  ihe  discretized  formulation.  It  is  also  valid  for  piecewise 
linear  strain-hardening. 

The  extrusion  thickness  has  been  observed  to  be  about  0.1  tttU.  (micron) 
1381 

by  "roT^'n  and  Stubbington  ,  so  that  the  distance  between  the  two  slices  is 
taker,  to  be  0.1  h)jA.  .  The  thickness  of  the  slices  P  and  Q  in  Figure  8  is  very 
small  and  is  assumed  to  be  0.01  mjA  .  The  Initial  stress  field  in  the  primary  slip 
system  at  the  free  surface  is  assumed  to  be  zero  everywhere  except  in  the  two 
neighboring  slices  P  and  Q,  in  which  the  initial  stress  varies  linearly  with  the 
depth  of  the  slice,  i.e., 

r  SO-  inP 

rjtx.o  =  /  lnQ 

k.  O  elsewhere 

where  the  linear  dimension,  d,  of  the  sliding  crystal  is  taken  to  be  25  mjA,  . 

To  reduce  the  lengthly  computer  calculation,  the  secondary  slip  is  allowed  to 
slice  only  In  the  region  between  two  slices  P  and  Q,  and  up  to  a  distance 
10  irijA.  from  the  free  surface.  200  grids  within  this  10  distance  were  taken 
for  this  secondary  siip  calculation.  The  sliding  of  the  secondary  slip  system  in 
the  inferior  region  between  P  and  Q,  and  beyond  10  mjA.  from  the  free  sur¬ 
face  is  taken  to  have  little  effect  on  the  plastic  strain  build-up  near  the  free 
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surface.  The  relief  of  the  resolved  shear  stress  due  to  plastic  strain  in  tha  slid  grid 

is  assumed  to  be  given  by  that  at  the  centroid  of  the  grid.  The  build-up  of  local 

plasric  strain  In  the  primary  and  secondary  slip  systems  is  calculated  for  an  applied 

alternating  stress  of  100  psi  along  x 2  -direction-  The  critical  shear  stress 

is  Taken  to  be  53.5  psi,  the  shear  modulus  ji.  to  be  3.85  x  10^  psi  and  the 

elastic  Poisson's  ratio  0”  to  be  0.3.  These  values  correspond  approximately  to 

'391 

those  for  pure  aluminum  1  . 

The  calculated  resolved  shear  stress  distributions  under  forward  loading 
at  tne  35th  cycle  with  and  without  considering  the  secondary  slip  are  shown  re¬ 
spectively  in  Figures  9  and  10.  It  is  interesting  to  note  that,  during  fatigue  loadings, 
the  difference  in  the  residual  resolved  shear  stresses  in  the  slices  P  and  Q  near 
the  free  surface  is  about  constant  for  the  case  with  only  primary  slip  system  sliding, 
wnlle  the  difference  is  steadily  increased  with  both  primary  and  £  econdary  slip 
systems  jointly  sliding. 

The  calculated  local  plastic  strain  in  the  primary  slip  system  in  the 
Slice  P  considering  sliding  In  the  secondary  siip  system  at  various  cycles  of  loading 
Is  shown  in  Figure  1 1 .  The  local  plastic  strain  reaches  approximately  3%  at  the 
vree  surface  grid  in  about  36  cycles.  The  local  plastic  strain  build-up  at  the 
.r.terior  end  of  the  thin  slice  is  much  slower  (only  about  1/3)  than  that  at  the 
free  surface  as  predicted  by  Lin  and  !to  .  .  Hence,  in  general,  fatigue  cracks 

occur  at  the  free  surface  before  any  interior  crack  has  a  chance  to  initiate. 

The  increase  of  maximum  local  plastic  strain  in  the  primary  slip  system 
at  tne  free  surface  with  and  without  the  sliding  of  the  secondary  slip  system  are 
suov/r.  in  Figure  1  2.  The  build-up  of  the  plastic  strain  in  the  primary  slip  system 
with  the  secondary  slip  sliding  Is  much  faster  than  that  with  only  primary  slip 
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system  siloing.  The  build-up  in  the  letter  cose  Is  approximately  linear  with  the 
number  of  cycles. 

V  CONCLUSION 

The  resolved  shear  stress  field  caused  by  a  uniform  distribution  of 
plastic  shear  strain  in  one  grid  in  a  semi-infinite  medium  was  obtained  in  closed 
form.  This  solution  in  closed  form  was  used  to  calculate  the  build-up  of  local 
plastic  strain  at  the  free  surface.  The  local  strain-hardening  in  the  slid  region  is 
assumed  to  be  zero.  The  activation  of  the  second  slip  system  is  allowed  during  the 
fatigue  cycling.  Calculated  results  show  that  sliding  of  the  primary  slip  system 
increases  rhe  residual  resolved  shear  stress  in  the  secondary  slip  system  monotonically 
at  a  relatively  rapid  rate  and  helps  the  secondary  slip  system  to  slide  at  the  early 
state  of  fatigue  cycling.  The  slip  in  the  secondary  slip  system  also  affect  greatly 
the  rate  of  the  local  plastic  strain  build-up  in  the  primary  slip  system.  The  build¬ 
up  of  iocal  plastic  strain  is  used  as  the  criterion  for  estimating  early  fatigue  damage. 
Hence,  the  sliding  of  the  secondary  slip  system  increases  greatly  the  rate  of  fatigue 
damage,  in  those  regions  where  two  slip  systems  are  active,  the  total  plastic 
strain  is  the  tensor  sum  of  the  slip  of  the  two  systems.  The  amounts  of  slip  in  these 
two  systems  vary  from  point  to  point.  The  plane  of  maximum  plastic  shear  strain 
caused  b)  slip  in  these  two  systems  varies  from  point  to  point.  It  seems  that  the 

variation  of  the  orientation  of  this  plane  may  cause  the  waviness  of  slip  lines 

,  .  [41,42] 

ooserved  on  specimen 

The  slip  distribution  calculated  by  the  present  analysis  satisfy  the 
conditions  of  compatibility  and  equilibrium,  as  well  as  the  dependency  of  slip  on 
rhe  resolved  shear  stress  throughout  the  metal  for  all  stages  of  alternating  load. 
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Part  C 


Preliminary  Evaluation  of  Plasma  Spray 
Metal  Coatings  as  a  Means  for  Prolonging  Fatigue  Life 

The  nucleation  of  a  fatigue  crack  generally  occurs  at  the  free  surface 
due  to  the  large  build-up  of  local  plastic  strain  resulting  from  slip  at  the  free  surface. 
It  was  experimentally  shown  that  successive  electropolishing  the  surface  of  a  specimen 
undergoing  fatigue  testing  will  considerably  increase  the  fatigue  life.  Honce,  it 
appears  that  the  fatigue  life  of  a  specimen  can  be  extended  if  its  surface  fatigue 
properties  are  improved. 

One  method  that  may  accomplish  this  is  to  give  the  specimen  surfa:<*  a 
coating  which  has  a  higher  fatigue  life  than  the  specimen  itself.  This  coating  must 
form  a  strong  molecular  bond  with  the  specimen.  It  was  thought  that  some  plasma 
spray  coatings  might  satisfy  these  requirements.  A  brief  literature  survey  was  made 
on  Ihe  effect  of  plasma  spray  coating  on  fatigue  life  of  base  material.  The  results 
of  this  survey  are  briefly  summarized  as  follows. 

Plasma  spray  coatings  of  Nlckel-Aluminide  and  Carbide  types  were  found 
to  adversely  affect  the  fatigue  properties  of  aluminum  alloy  2024-T4  and  ultra-high 
strength  steel,  AISI  type  H-ll  The  plasma  spray  process  or  preparation  procedure 

has  little  effect  on  the  high  cycle  endurance  limit  of  aluminum  alloy  2024-T4.  This 
agrees  with  previous  work  done  by  Whittaker  Remmelts  and  Wolff 
However,  there  is  some  low  cycle  fatigue  life  reduction.  The  high  cycle  endurance 
limit  of  high  strength  H-l  1  steel  heat-treated  to  280  ksi  tensile  strength  was  found 
to  reduce  by  15  to  25  percent  by  tho  application  of  plasma  sprayed  Nlckel-Aluminide 
and/or  Carbide  coatings,  with  corresponding  reduction  in  low  cycle  fattgue  life. 

The  same  results  were  found  In  the  plasma  coated  titanium  6AL-4V  and  AISI  9310 
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The  above  plasma  spray  coatings  do  not  have  beneficial  effects  on  fatigue 
However,  the  fatigue  lives  of  the  coating  materials  under  cyclic  strain  are  not 
known.  The  reduction  of  fotlnue  life  of  the  base  marerial  due  to  the  coating  may 
be  caused  by  the  low  fatigue  life  of  the  coating*  from  Coffin- Meson's  law  of  . 
fatigue,  under  low  cycle  fatigue  loadings,  ductile  materials  nave  longer  fatigue 
life  under  constant  cyclic  strain  than  the  less  ductile  ones.  A  coating  of  ductile 
material  well  bonded  to  a  less  ductile  base  should  increase  the  fatigue. life  of  the 
base  material .  However,  such  a  metallurgical  process  giving  perfect  bonding  pi 
this  coating  to  such  o  base  material  has  not  been  found  in  the  literature.  For  the 
Improvement  of  fatigue  properties  such  o  process  needs  to  be  developed. 
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APPENDIX 


Let  the  plastic  strain  gradient  In  slid  region  S  be  written  as: 

n..  m  .  -  _  ..  .. 

•£57"—"  QO*-*i  v  for  xeS  o*  x-»  iT  (A-l) 

where  p  Is  a  real  number.  For  P-"l,e"(x,m)  would  give  a  logarithmic 
singularity  and  becomes  Infinite  as  x  approaches  x.  Consider  plastic  strain  as 
caused  by  the  displacement  of  dislocations.  Let  the  lattice  spacing  be  denoted 
by  "A1*  which  is  a  finite  distance.  The  change  of  plastic  strain  across  "A"  Is 
always  finite.  Hence  Is  bounded  and  p  In  Equation  (Arl)  Is  greater 

than  zero.  The  stress  Green's  function  Z^  (X,Z)  at  x  due  to  a  unit  force  noting 
at  x  along  x^  direction  has  been  shown  to  be  related  to  the  stress  function 
(Equation  3),  It  is  easily  observed  through  the  expressions  for  the  stress  function 
that 

Zy  mm  0  <  ix-xf)  for  X4  S  as  x—5?  (A-2) 

Let  S,,(x)  be  a  sphere  on  ( x j,  Xj)  plane,  centered  at  x  with  radius  £  .  It  can 
be  shown  from  Equations  (A-1 )  and  (A-2)  that  for  x  0  (x) 

|  as  e-~o  (Ar3) 

Sjiu  k,  t  —  x 

This  integral  approaches  zero  for  p  >-/.  Since  p> 0 ,  this  integral  vanishes  as 
.  Even  though  the  stress  Green's  function  Zn  (**% )  Is  singular  at 

x  -  x,  the  Integral  J  Vy  <KtJ  converges  uniformly  In  the 

slid  region  S, 
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